Magnetism of Fe and Ni from rotationally invariant 
Hirsch-Fye quantum Monte Carlo calculations 
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We present a computational scheme for the investigation of complex materials with correlated elec- 
trons which allows one to take into account the spin rotational symmetry of the exchange Coulomb 
interaction. It combines ab initio band structure methods with dynamical mean-field theory and 
is implemented with the Hirsch-Fye quantum Monte Carlo algorithm in which the spin rotational 
invariance of Hund's exchange is restored by averaging over all possible directions of spin quantiza- 
tion axis. We employ this technique to perform benchmark calculations for the three-band Hubbard 
model on the infinite-dimensional Bethe lattice. Our results agree quantitatively well with those 
obtained using the continuous-time quantum Monte Carlo method. The proposed approach is em- 
ployed to compute the electronic and magnetic properties of paramagnetic a iron and nickel. The 
obtained Curie temperatures agree well with experiment. Our results indicate that the magnetic 
transition temperature is significantly overestimated by using the density-density type of Coulomb 
interaction. 

PACS numbers: 71.15.Mb, 71.20.Be, 71.27.+a 



I. INTRODUCTION 

The theoretical description of the electronic proper- 
ties of transition metal compounds with partially filled d 
and /-shells and strong Coulomb interaction between the 
electrons remains a challenging, fundamental problem in 
condensed matter physics The interplay between elec- 
tronic and lattice degrees of freedom in such materials 
results in their diverse physical properties and rich phase 
diagrams making these compounds particularly attrac- 
tive for technological applications^ Moreover, orbital de- 
generacy is an important and often inevitable cause of 
this complexity. Together with the Hund's exchange in- 
teraction, it has important implications for the electronic 
and magnetic properties of correlated materials, leading 
to formation of local moments and complicated multiplct 
structures. 

The electronic properties of correlated materials can 
be understood by employing the so-called LDA+DMFT 
approachj^ a combination of ab initio local density ap- 
proximation (LDA) of the density functional theory and 
dynamical mean- field theory (DMFT). Nowadays, the 
LDA+DMFT technique has become a state-of-the-art 
method for realistic description of correlated electron ma- 
terials from first principles. This approach provides a 
systematic many-body treatment of the effect of local 
electronic correlations by taking into account temporal 
fluctuations while spatial fluctuations arc neglected. Ap- 
plications of LDA+DMFT for correlated electron com- 
pounds such as transition metals and their oxides have 
provided important insights into our understanding of 
the electronic and magnetic properties of these materi- 
als. In particular, by employing the LDA+DMFT tech- 
nique it has become possible to obtain a good quantita- 



tive description of localized as well as delocalized electron 
states. In addition, the approach allows one to determine 
the electronic and magnetic properties of correlated com- 
pounds in both paramagnetic and magnetically ordered 
states. 

Nevertheless, there are two important limitations 
of conventional implementations of the LDA+DMFT 
method. The first originates from the single-site (local) 
nature of DMFT. In particular, the key assumption of the 
theory is the limit of infinite spatial dimension, which 
allows one to perform an exact mapping of a complex 
lattice model (such as the Hubbard model) to a quan- 
tum impurity with an energy-dependent external bath, 
resulting in fc-independent self-energy. However, in some 
cases the non-local spatial correlations can be essential to 
provide correct description of the properties of correlated 
materials^ For instance, the standard LDA+DMFT cal- 
culations are not able to capture the reduction of mag- 
netic transition temperature due to long- wavelength spin 
waves. To resolve this problem there have been recently 
proposed several methods^ which we leave beyond the 
scope of our paper. 

Since correlated materials often have several bands at 
the Fermi level, it requires specific treatment of the lo- 
cal Coulomb interaction, which in a cubic environment 
consists of the intra- and intcrorbital Coulomb interac- 
tions U and U', Hund's exchange J, and the pair-hopping 
coupling J'. These interactions obey spin and orbital ro- 
tational symmetry, thereby U = U' + J + J' ensures the 
rotational invariance in the orbital space and J = J' can 
be assumed whenever the spin-orbital coupling is neg- 
ligible. Unfortunately, it is difficult to handle all these 
multiband interactions including Hund's exchange cou- 
pling and the pair hopping term with the quantum Monte 
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Carlo (QMC) method. In particular, a straightforward 
implementation leads to a severe sign problem making 
such simulations unfeasible. Therefore, at present, the 
most material-specific calculations employ the approxi- 
mate form of the Coulomb repulsion restricted to the 
Ising-type exchange interaction. These calculations often 
provide a good quantitative description of the electronic, 
magnetic, and structural properties of correlated materi- 
als as a function of the reduced temperature T/T m , where 
T m is the calculated temperature of magnetic ordering^ 
However, the correct symmetry of the exchange interac- 
tion turns out to be essential for quantitative description 
of the electronic and magnetic properties of correlated 
systems.— ~— 

This problem can be overcome by using some quantum 
impurity solvers such as numerical renormalization group 
(NRG)r^ exact diagonalization (ED)^ continuous-time 
quantum Monte Carlo (CT-QMC)^ and others,!^ 
which allow one to treat the Coulomb interaction in its 
general form with preserved spin rotational symmetry. 
These calculations performed for the two- and three-band 
Hubbard models on the infinite-dimensional Beth o 10 i 12 
and hypercubioii lattices show substantial overestima- 
tion of the magnetic transition temperature for the ap- 
proximate Ising-type form of the exchange Coulomb in- 
teraction with respect to the rotationally invariant one. 
In accord with this, recent LDA+DMFT calculations of 
correlated compounds also indicate that the magnetic 
transition temperatures appear to be significantly overes- 
timated by using the density-density type of Coulomb in- 
teraction.—^^ However, applications of these techniques 
so far have been mostly limited to simple model systems 
and only a few realistic calculations for the 3d compounds 
have been recently presented^ This is mostly because of 
high computational costs (exponential with the number 
of orbitals) of these methods implemented with the full 
rotationally invariant Coulomb interaction which makes 
such calculations for the 3d and 4/ materials extremely 
expensive. Obviously, the LDA+DMFT investigations of 
correlated materials with the Coulomb interaction in its 
general form with preserved spin rotational symmetry re- 
main problematic and pose a great theoretical challenge. 

In this paper, we present a novel implementation of 
the LDA+DMFT approach which allows us to take into 
account rotational symmetry of the exchange Coulomb 
interaction. This approach is formulated in terms of the 
Hirsch-Fye quantum Monte Carlo algorithm^ where the 
spin rotational invariance of Hund's exchange is restored 
by averaging over all possible directions of spin quantiza- 
tion axis. It provides a robust and computationally effi- 
cient method which allows us to simulate the five-orbital 
systems at high temperatures. Using this technique we 
perform benchmark calculations for the three-band Hub- 
bard model on the infinite-dimensional Bcthe lattice. In 
addition, we employ the proposed approach to calculate 
the electronic and magnetic properties of paramagnetic a 
iron and nickel. To outline the importance of rotational 
symmetry of the exchange Coulomb interaction we com- 



pare our results with those obtained by using the density- 
density approximation. 

This paper is organized as follows. In Sec.|H]we present 
a detailed formulation of the proposed approach which 
allows one to treat rotational invariance of the exchange 
interaction. In Sec. lIIII we employ this technique to com- 
pute the electronic and magnetic properties of the three- 
band model on the Bethe lattice, paramagnetic a iron, 
and nickel. The obtained results are compared with those 
of previous calculations and experimental data. Finally, 
conclusions are presented in Sec. IIVI 



II. METHOD 

The multiband Hamiltonian with full rotationally in- 
variant on-site Coulomb interaction^ can be written in 
the following form 

H = U ^ "mt"rr4 (!) 
m 

+ - ^2 i( u ~ 2J)h ma h m icr + (U - 3J)n ma n m 'a 

mm/a 
m^m 

where c+ CT (c mcr ) denotes the creation (annihilation) op- 
erator of an electron with spin a (=t, i) at orbital m, 
hma = ^mcfimai U is the screened Coulomb interaction 
parameter, J is the Hund's exchange coupling. The 
first three terms in Hamiltonian ([T]) correspond to the 
density-density part of Coulomb interaction and contain 
the exchange interaction in the Ising-type form. The re- 
maining part consists of spin-flip (4th) and pair hopping 
(5th) terms. Using the z-projection of the spin operator, 
Sj n = (n m ^ — n m j_)/2, and the orbital occupancy opera- 
tor, N rn = n rn ^ + fimi, the density-density part can be 
rewritten as 

Hdd = U^^ hmfhmi (2) 

m 

+ \Y1 {UN m N m '~2JS^S z m ,}, 

mm' 
m^m' 

where U = U — 5 J/2 is the average value of the Coulomb 
interaction. The spin-flip term in Eq. (JT]) can be ex- 
pressed via operators = (cj^c^ + & m ]C m ^,)/2 and 

&m = — *(^mt*W — c m ±c m ^)/2 as 

^ ] ^mtr^ma^m/W^mla ~ ^(^m^m' ^m^m')' (*0 
a 

The pair hopping term acts only on high energy states 
with two electrons on the same orbital and thereby can 
be neglected. However, taking into account the spin- 
flip term in Eq. ([T]) is crucial for the correct description 
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of spin dynamics. Therefore, the Coulomb interaction 
Hamiltonian can be written as 



(4) 



1 



2 ^ 

mmt 



{UN m N m , -2JS m S m .}. 



The Hamiltonian (4) with exchange term taken as a 

vector product, JS m S m > , is invariant with respect to the 
spin quantization axis rotations while the density-density 
counterpart with the Ising-type exchange term, JS^S^,, 
is not. 

To restore the spin rotational symmetry of Hamilto- 
nian we here employ the method originally proposed 
by Hubbard^i The Coulomb interaction in Ref. [2l| was 
considered in the following form 

UN t N i = -UN 2 -US 2 = ^UN 2 -U(eS) 2 , (5) 



where N a = J2, 



N = N t + N u S z = J2 m S? n ,Sk 



the total spin of atom, e is an arbitrary unit vector that 
can be interpreted as a quantization axis. In order to 
restore the spin rotational symmetry in Eq. (|5|) , averaging 
over all directions of e was introduced. 

In this paper, we propose to restore the spin rotational 
invariance in Eq. ([2]) by averaging in the partition func- 
tion calculation over all possible directions of spin quan- 
tization axis determined by the angles <fi and 9: 



Z oc / d9d(j)Tt[e 



-P(Ho + U 1 (e,<p)H dd U(9^)h 



(0) 



Here, Hq is a non-interacting part of impurity Hamilto- 
nian, U ac ,' (9, </>) stands for a transformation matrix for 
spin variables and reads as 

U(9 *,)-( cos(#/2)e^ 2 sin(0/2) e ^/ 2 \ 
U{U ^> ~ \-sin{9/2)e^/ 2 cos(# /2)e~^ 2 J' [t) 

According to Eq. [6j result does not depend on the choice 
of quantization axis direction and hence is rotationally 
invariant. 

In the Hirsch-Fyc quantum Monte Carlo algorithm^ 
dynamical potential matrix V ma {T) is diagonal in spin 
and orbital indices and can be expressed as 

1, [L <V 
-1, /i > V ' 

where [i and v are combined spin-orbital indices, 
s M „ is an Ising-likc variable taking the values ±1, 

•Vi, = arcosh(e ArC/ ^"/ 2 ). To average over all possible di- 
rections of spin quantization axis in the partition function 
(Eq. [6J , one should consider the following potential 



V l {r) = U\e,<j>)V{T)U{6,4>), 



(9) 



which is non-diagonal in spin index. Therefore, in QMC 
simulation Markov's chain should be defined not only 
by stochastic change of fluctuating fields but also by the 
stochastic change of angles 9 and 0, corresponding to the 
quantization axis direction. 

The physical meaning of the proposed rotationally in- 
variant algorithm can be expressed as averaging over all 
possible directions of fluctuating spin polarization in the 
3d shell, in contrast to polarization along the z-axis only 
in density-density Hamiltonian ([2]). Thereby this tech- 
nique allows one to take into account not only the lon- 
gitudinal spin fluctuations but also the transverse ones, 
reproducing the full spin dynamics. 



III. RESULTS AND DISCUSSION 

In this section, as a proof of concept we first present 
results of our model calculations in comparison with pre- 
vious studies. In particular, we employ the proposed ap- 
proach to compute the three-band Hubbard model on the 
infinite-dimensional Bethe lattice. We benchmark these 
calculations with the previously published results of the 
CT-QMC computations^ To study the role of symmetry 
of the exchange interaction we compare our results with 
those obtained by employing the density-density approx- 
imate form of the local Coulomb interaction. 

To proceed further, we investigate the electronic and 
magnetic properties of paramagnetic iron and nickel. We 
outline importance of the correct spin and orbital ro- 
tational symmetry of the exchange interaction on the 
magnetic properties of these materials. To compute the 
magnetic properties of correlated materials we evaluate 
uniform magnetic susceptibility x(T) = dM(T)/dH z by 
measuring the induced magnetization in a small external 
magnetic field. Here, M(T) = X)m( n ™t — n mi) is the 
magnetization induced by the external magnetic field H z 
applied along the z axis (these calculations include the 
polarization of the impurity Weiss field due to the exter- 
nal magnetic field H z ). In our calculations, we adopted 
a few magnetic fields in the range from 0.01 to 0.04 eV 
to ensure linearity of response. 

The calculated inverse uniform magnetic susceptibility 
exhibits a linear behavior at high temperatures which 
indicates the presence of local magnetic moments. The 
temperature dependence of x(T) obeys the Curie- Weiss 
law x(T) = C/(T — Tc), where C is a material-specific 
constant and Tq is the Curie temperature. The result 
of the least-square fit to the Curie- Weiss law is shown 
by the straight line. Therefore, the Curie temperature 
can be determined as an intersection of this line with the 
temperature axis. 



A. Three-band model 

In recent years the properties of the two- and three- 
orbital Hubbard models have been extensively inves- 
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tigated by using dynamical mean-field approach^ -11 ' 22 
Here we would like to refer to Ref. [l^ where the three- 
orbital Hubbard model on the infinite-dimensional Bet he 
lattice was studied by means of the continuous-time 
quantum Monte Carlo method with full rotationally in- 
variant Coulomb interaction. In particular, it was es- 
tablished that the effect of spin-flip interactions consid- 
erably depends on a band filling. It was shown that the 
Mott-Hubbard physics dominates due to the strong effec- 
tive Coulomb interaction U in the particle-hole symmet- 
ric case, while away from half-filling formation of local 
magnetic moments is more plausible due to the Hund's 
exchange. For the benchmark purposes we here repro- 
duce these calculations by using the rotationally invari- 
ant HF-QMC method. We have chosen the same set of 
parameters as in Ref. Oil Namely, we consider the three- 
band model with two-electron occupancy. We adopt 
U = 8t and J = 1.2t for the Coulomb interaction and 
Hund's exchange, where t is a half of the non-interacting 
bandwidth. In Fig. [1] we present our results for the in- 
verse uniform magnetic susceptibility. The calculated 
Curie temperatures are 0.29i and 0.1 5i for the conven- 
tional HF-QMC and rotationally invariant methods, re- 
spectively. These findings are in good quantitative agree- 
ment with results of the CT-QMC calculations, which 
give Tc ~ 0.27i and 0.14i for the density-density and 
full rotationally invariant exchange, respectively. This 
quantitative agreement clearly demonstrates the validity 
of our method for accurate description of the spin dy- 
namics. We note that transverse spin fluctuations can be 
regarded as an important source of softening of magnetic 
response. It is expected that in the five-band case this 
effect should be even more profound. 



B. a iron 

Elemental iron is one of the oldest and experimentally 
best studied itinerant ferromagnets. Various properties 
of iron can be understood on the basis of band-structure 
calculations.— In particular, these calculations provide 
a good description of the low-temperature ferromagnetic 
phase of Fe. However, applications of conventional band- 
structure techniques to describe the properties of param- 
agnetic iron, in particular, close to the a — 7 phase tran- 
sition, do not lead to satisfactory results. This is mainly 
due to the presence of local magnetic moments above 
the Curie temperature which are important for quanti- 
tative description of paramagnetic state. In this respect, 
the LDA+DMFT approach provides the best formalism 
which allows one to unify the localized and itinerant elec- 
tron behavior in metallic magnets. 

Applications of LDA+DMFT have shown to provide 
a good quantitative description of the electronic, mag- 
netic, and structural properties of iron . 7 ' 8 ' 24 " — However, 
an agreement was achieved only in terms of the reduced 
temperature T/Tc, while the calculated Curie temper- 
ature Tc was found to be almost twice larger than the 
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FIG. 1. (Color online) Temperature dependence of the in- 
verse uniform magnetic susceptibility obtained by the DMFT 
for the three-band infinite-dimensional Bethe lattice. The pa- 
rameters are chosen to be the same as in Fig. 3 of Ref. fig . 
The straight lines depict the least-squares fit to the Curie- 
Weiss law. The obtained Curie temperature values are 0.29* 
and 0.15* for HF-QMC and rotationally invariant HF-QMC, 
respectively. Our findings are in good quantitative agreement 
with results of the CT-QMC calculations presented in Ref. Hcl. 



experimental value of 1043 K (Ref. I29I ). Recently the 
properties of iron have been investigated by means of J- 
QMC method, which uses the static approximation for 
the charge degrees of freedom and treats the exchange 
Coulomb interaction in the rotationally invariant form^, 
These calculations indicate that substantial part of the 
Curie temperature overestimation comes from the ap- 
proximate (density-density) treatment of the exchange 
Coulomb interaction. 

We now calculate the electronic structure and mag- 
netic properties of paramagnetic bec iron by employing 
the LDA+DMFT implemented with the rotationally in- 
variant HF-QMC method. We first calculate the non- 
magnetic LDA electronic structure of a iron using the 
tight-binding linear muffin-tin orbital (TB-LMTO) ap- 
proach^ For these calculations we adopt lattice constant 
a = 2.866 A. We construct an effective low-energy Hamil- 
tonian in the basis of Fe spd Wannier orbitals using the 
Ath-order muffin-tin (AMTO) method^ In accord with 
the previous estimations,— ' 24 ' 25 ' 32 we take U — 2.3 eV 
and J = 0.9 eV for the screened Coulomb interaction 
and Hund's exchange, respectively. 

In Fig. [5] we present the partial density of states and 
the corresponding imaginary parts of the self-energies ob- 
tained by the rotationally invariant HF-QMC method at 
T = 1160 K in comparison with the non- magnetic LDA 
results. Our calculations reproduce the splitting in the 
density of states of the e g orbitals near the Fermi level, 
which is absent in the LDA calculations. This result 
agrees well with the previous calculation s 13 ' 24 ' 25 as well 
as with the experimental data and is one of the char- 
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FIG. 2. (Color online) Partial t2 g (black) and e g (red) den- 
sities of states for paramagnetic a iron as obtained by the 
rotationally invariant HF-QMC within LDA+DMFT in com- 
parison with the non-magnetic LDA results. The Fermi level 
is indicated by the vertical (gray) line at zero energy. Inset: 
imaginary parts of the obtained self-energies. 



acteristic features of a iron. The calculated self-energy 
for the tig states exhibits a Fermi-liquid-likc behavior, 
whereas the e g self-energy diverges at low frequencies. 
Our calculations with J = recover the Fermi-liquid-likc 
behavior for the e g states resulting in the suppression of 
the splitting. This indicates that the splitting can be 
attributed to the exchange Coulomb interaction^ We 
found no evidence for the Hubbard subbands formation 
in the calculated quasiparticle spectrum. Thereby we 
conclude that in a iron the correlation effects are mainly 
effected by the strength of the Hund's coupling J rather 
than by the Coulomb interaction U which allow us to re- 
fer a iron as a Hund's metal. We note that importance 
of the Hund's exchange in multiorbital systems has been 
recently studied in Ref. |H. 

The temperature dependence of the inverse uniform 
magnetic susceptibility calculated by the LDA+DMFT 
shows a linear behavior at high temperatures (Fig. [3]) in 
accord with the Curie- Weiss law. It is clearly seen that 
the HF-QMC limited to the Ising-type exchange interac- 
tion substantially overestimates the Curie temperature 
value and yields Xc ~ 2050 K. The rotationally invari- 
ant HF-QMC method gives T c ~ 1260 K, which is in 
good agreement with the experimental value of 1043 K 
(Ref. |29I ). The calculated values of effective local mo- 
ments extracted from the uniform magnetic susceptibil- 
ity are /x cff ~ 2.87 fi B and 2.61 fi B for the HF-QMC with 
density-density and rotationally invariant algorithms, re- 
spectively. Our estimates are in reasonable agreement 
with the experimental value of 3.13 /xb (Ref. |29( ). which 
appears to be underestimated by both methods. 

We note that preserving the spin rotational symme- 
try turns out to be crucial for quantitative description of 
magnetic properties of correlated materials. In particu- 
lar, for a iron this leads to substantial improvement of the 



magnetic transition temperature value. The remaining 
discrepancy can be attributed to non-local effects which 
are not present in conventional DMFT. 



C. Nickel 

Elemental nickel is another example of itinerant elec- 
tron ferromagnets which together with iron serves as 
a benchmark material for electronic structure methods. 
Various low temperatures properties of nickel can be 
understood by employing standard band-structure ap- 
proaches^ Nevertheless, these techniques are generally 
fail to reproduce many characteristic features of nickel, 
such as an existence of satellite structure^ at —6 eV, 3d 
electron bandwidth,— and the value of exchange split- 
ting^ Applications of LDA+DMFT to study the elec- 
tronic and magnetic properties of nickel have given a 
good quantitative description of many of these phenom- 
enai 7 ' 25 ' 27 ' 28 ' 37 ' 38 The calculated magnetic properties of 
nickel are shown to be in good agreement with experi- 
ment. In contrast to iron, the overestimation of the mag- 
netic transition temperature for nickel by LDA+DMFT 
with the Ising-type exchange interaction is not so signifi- 
canti 7 -^ However, as demonstrated below, preserving the 
spin rotational symmetry leads to the underestimation of 
the Curie temperature. 

We now compute the electronic structure and mag- 
netic properties of paramagnetic nickel by employing the 
LDA+DMFT implemented with rotationally invariant 
HF-QMC method. To obtain the non-magnetic LDA 
electronic structure of fee nickel we employ the tight- 
binding linear muffin-tin orbital approach. The calcula- 
tions were performed for the lattice constant a = 3.524 A. 
Using these results we construct an effective low-energy 
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FIG. 3. (Color online) Temperature dependence of the inverse 
uniform magnetic susceptibility for 01 iron as obtained by the 
LDA+DMFT. The straight lines depict the least-squares fit to 
the Curie- Weiss law. The experimental value Tq XP = 1043 K 
is denoted by the (black) arrow. The experimental value of the 
effective local magnetic moment is fj!^ p = 3.13/iB (Ref. [2^). 
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Hamiltonian in the basis of Ni spd Wannicr orbitals by 
employing the AMTO method. In accord with the previ- 
ous estimations,— we take U = 2.3 eV and J = 1.0 eV for 
the screened Coulomb interaction and Hund's exchange, 
respectively. 

In Fig. [4] we present the partial densities of states and 
the imaginary parts of the self-energies obtained by the 
rotationally invariant algorithm at T = 1160 K in com- 
parison with the non-magnetic LDA results. The inclu- 
sion of electronic correlations results in a small reduc- 
tion of the bandwidth of nickel with respect to the non- 
magnetic LDA result. In addition, a satellite-like struc- 
ture emerges at about —5.5 eV. In contrast to iron, the 
obtained self-energies for both the tig and e g orbitals ex- 
hibit the Fermi-liquid-like behavior. This can also be seen 
from the calculated amplitudes of the effective damping 
3E(0) which are -0.03 eV and -0.02 eV for the t 2g and 
e g states of nickel, respectively, while Sffi(0) ~ —0.24 eV 
for the t2 g states of a iron. This indicates more coherent 
nature of the electronic properties of nickel in compari- 
son with iron. To quantify this qualitative differenc e 7 ! 39 
we calculate spin-spin correlation functions for iron and 
nickel. In Fig. [5] we present the impurity spin-spin cor- 
relation functions on the real and imaginary energy axes 
calculated for a iron and nickel at T = 2.5 Tq (Tq refers 
here to the corresponding calculated value of the Curie 
temperature). The height of peak on the real energy axis 
can be interpreted as a value of the local magnetic mo- 
ment. The pronounced peak for iron indicates the pres- 
ence of well localized magnetic moments above Tc, while 
magnetism of nickel is more itinerant. 

In Fig. |5] we show the temperature dependence of 
the inverse uniform magnetic susceptibility calculated by 
LDA+DMFT. From these data we estimate the values 
of the Curie temperature. By employing the density- 
density approximation we find Tq ~ 840 K, whereas the 
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FIG. 4. (Color online) Partial densities of states for paramag- 
netic Ni obtained by rotationally invariant HF-QMC within 
LDA+DMFT in comparison with the LDA ones. The Fermi 
level is indicated by the vertical (gray) line at zero energy. 
Inset: imaginary parts of the obtained self-energies 
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FIG. 5. (Color online) Spin-spin correlation functions on the 
real and imaginary (inset) energy axes for a iron and nickel 
calculated by LDA+DMFT at T = 2.5 T c . 



inclusion of the spin rotational symmetry leads to al- 
most twice smaller value of about 400 K. Both results 
are in reasonable agreement with the experimental value 
of 631 K. To proceed further, we compute the effective lo- 
cal magnetic moments by the HF-QMC and rotationally 
invariant methods which give /i c g ~ 1.55 /Ub and 1.49 /xb, 
respectively. These findings are in good agreement with 
the experimental value of 1.62 [is (Ref. [291 ) . The fact 
that in a iron the LDA+DMFT with proper spin rota- 
tional symmetry results in an overestimated value of Tc, 
while it appears to be underestimated in nickel, can be 
dealt with non-local effects which arc neglected in DMFT 
or with the more itinerant nature of magnetism in nickel 
than in a iron. 
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FIG. 6. (Color online) Temperature dependence of the in- 
verse uniform magnetic susceptibility for Ni as obtained by 
LDA+DMFT. The straight lines depict the least-squares fit 
to the Curie- Weiss law. The experimental value T£ xp = 631 K 
is denoted by the (black) arrow. The experimental value of the 
effective local magnetic moment is /i^ p = 1.62 /j,b (Ref.[2S|). 
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IV. CONCLUSIONS 



In conclusion, we have presented a novel implemen- 
tation of the LDA+DMFT approach which allows one 
to take into account the spin rotational symmetry of 
the exchange Coulomb interaction. The computational 
scheme is based on extension of the Hirsch-Fye quantum 
Monte Carlo algorithm in which the spin rotational in- 
variancc of Hund's exchange is restored by averaging over 
all possible directions of spin quantization axis. The pro- 
posed approach provides a robust and computationally 
efficient method which allows us to compute high tem- 
perature electronic properties of the five-orbital systems. 
We have used this approach to perform benchmark calcu- 
lations for the three-band Hubbard model on the infinite- 
dimensional Bcthc lattice. Our results agree quantita- 
tively well with those obtained using the continuous-time 
quantum Monte Carlo technique. The proposed method 
is employed to compute the electronic and magnetic prop- 
erties of paramagnetic a iron and nickel. The obtained 
Curie temperatures agree well with experiment. Our re- 



sults indicate that the density-density approximation for 
the exchange Coulomb interaction leads to substantial 
overestimation of the magnetic transition temperature. 
The remaining discrepancy in the calculated Curie tem- 
peratures can be attributed to non-local effects which 
are not taken into account in standard implementations 
of the LDA+DMFT. 
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